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IfJinECHEiai  NKIVUKK  MXEDL  AND  HEURISTIC  FOR  IHE 
OCMBAT  SERVICE  SUPPORT  SUPPLY  SYSTEM  (CS4) 

INTRODUCnCN 


A  series  of  brainstoniilng  sessions  was  held  b/  the  Quaxtennaster 
SciiOQl  (QE)  through  the  Joint  Technical  Staff  (JTS)  at  Natick  Research, 
Developnent  and  Bigineering  Center.  The  involvement  of  the  Fcxxi 
Ehginei^ing  Directorate  (FED)  in  the  brainstoming  sessions  led  to  serious 
speculation  that  the  develoEsient  of  a  sultiechelon  inventory  model  of  a 
oosbat  service  support  supply  system  (CS4)  may  have  the  potential  to 
enhanoe  the  logistical  capabilities  in  any  given  theater  of  operations. 

In  any  theater  of  operations,  the  possibility  of  a  military  conflict 
is  always  present.  Not  only  oust  the  siqpply/disbribution  system  of  a 
theater  be  able  to  react  effectively  to  the  needls  of  the  troops,  it  must 
also  have  the  OE^)ability  to  sustedn  an  adeqjate  flow  of  sipplies  to  the 
field  troops,  in  a  coordinated  maimer,  in  the  event  of  a  prolonged  and 
intense  conflict.  Since  the  infomation  obtained!  by  intelligence 
gathering  is  not  always  accurate,  ocnplete,  or  current  enough  under 
constantly  changing  circumstances,  an  adequate  degree  of  prepareciness  and 
sustainability  is  essential  for  the  security  of  a  theater.  Without 
adequate  logii^cal  planning  and  siipport,  a  sudden  and  intense  conflict 
can  prevent  needed  subsistence  from  reaching  field  troops.  Particularly 
critioal  is  the  need  for  field  troop  sipply  support  at  the  lower  echelons 
txxm  the  \:pper  echelons  of  the  theater.  Not  only  must  the 
scpply/distribution  system  be  able  to  react  effectively  to  the  needs  of 
the  trxxps,  it  must  also  have  the  capability  to  sustadn  an  adequate  flow 
of  sqpplies  to  the  field  troops,  in  a  coordinated  manner,  in  the  event  of 
a  prolonged  and  intense  conflict. 

Supply/distribution  policies  that  are  currently  in  effect  in  theaters 
tend  to  be  static;  that  is,  they  do  not  account  for  the  stochastic  and 
dynamic  demands  that  can  adversely  influence  a  sipply/distribution  system 
in  times  of  military  conflict.  Thus,  under  the  current  policies,  the 
likelihood  of  problems  in  the  logistics  ^stem  that  are  severe  enou^  to 
ocBpaxmise  the  military  effectivaiess  of  a  theater  is  significant. 
Logistics  problems  that  are  of  paramount  concern  relate  to  the  possibility 
of  stookouts  oocurring  at  any  echelon  in  the  theater  as  the  result  of 
sudden  surges  in  demand;  these  are  e^secially  critical  at  the  lowest 
echelons  where  the  need  for  supplies  is  typically  most  highly 
cxnoentrated. 

There  are  many  analytical  methods  available  to  cope  with  logistics 
problems  that  are  related  to  stochastic  and  tynamlc  behaviors,  such  as 
desand.  Cperaticns  research  hzis  many  basic  methodologies;  however,  the 
method  or  oosbination  of  methods  applied  is  often  influenoed  by  judgment, 
perspective,  and  goals  of  the  anedyst.  Possibilities  that  exis^  within  a 
broad  range  of  tedmiques  include  inventory  analysis,  mathematicsd 
progranming,  just-in-time  inventory  methods,  and  network  analysis. 
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Ihe  objective  of  this  effort  is  to  identify  the  iisportant  aspects  of 
logisticzd  scpply/distributicn  probleois  vith  respect  to  military  theaters 
of  operation  develop  an  anedytical  e^proach  that  can  be  used  to  ensure 
that  troop  danands  for  perishable  rations  are  met.  The  goal  is  to  reduce 
inventory  levels  so  that  the  rations  are  rotated  more  guicikly  and  to 
reduce  the  use  of  tran^xartation  resources,  i.e. ,  v^iiides,  fuel, 
roadways,  railways,  etc. 

A  nebwoiic  analysis  approach  is  used  to  cut  inventory  levels  and 
transportation  atctivities  throughout  the  theater.  The  method  of  cfxsloe  is 
a  (Synonlc  flows  eopansion  of  a  minimum-cost  flow  method.  The 
out-of-Xilter  (QFA)  algorithm  is  used  primarily  because  lower  bounds  for 
the  arc  (route)  flows  not  only  ensure  geographical  oocundination  in  the 
logistics  netwozic  but  adso  enforoe  tesporal  coordination  as  well. 

Solution  methods  to  exogenous  source  limitations  and  pre^nsltioning  can 
also  be  aooounted  for  by  use  of  the  method  of  dynamic  flows. 

Since  a  major  aspect  of  the  modelling  is  the  existence  of  different 
types  of  rations,  the  flows  are  multiooDnodiiy.  A  heuristic  is  developed 
to  deed  with  the  nultiaoranodity  flows  version  of  the  problem. 
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THEATER  RATION  DISTKEBOTICN  SYSTEM 


This  section  of  the  report  outlines  the  adjects  of  the  ration 
sifply/distribution  system  that  lead  to  the  development  of  the  logistics 
modkL.  The  inportant  adjects  of  the  system  2u»  its  operating  policies, 
its  structure,  and  its  problematic  issues.  The  implication  of  the  mnrtei 
is  also  described.  Information  about  the  ration  distribution  systems  in 
military  theaters  was  drawn  from  ruaerous  discussions  with  eoperienoed 
personnel  and  ftxn  field  manuals. 

Qjeratina  Policies 

current  ration  distribution  policies  inclxide  a  "push”  and  a  "pull" 
ocnponent.  The  push  oonponent  projects  future  demands  and  moves  supplies 
aheiKl  in  an  effort  to  satisfy  the  projected  demands.  The  push  ocnponent 
moves  sipplies  (from  both  CONUS  and  OOCNUS)  to  facility  distribution 
locations,  currently  the  supplies  are  moved  by  long'-range  transportation 
to  any  given  theater;  this  is  usually  done  by  sea  in  30-day  pedletized 
lots.  The  lots  are  then  moved  by  thi^ter  transportation  tystems,  which 
are  tmuiaiy  truck  anVor  rail,  to  the  distribution  facilities  for  future 
regoirements.  The  pull  oonponent  of  the  distribution  process  reacts  to 
iaraediate  demands  that  occur  at  the  unit  level,  that  is,  primarily  at  or 
near  the  front  line  of  troops  (FLOT) .  Within  the  oonteset  of  this  project, 
the  term  "unit"  refers  generically  to  military  personnel  at  the  brigade 
level  or  lower  (see  Figure  1) .  Ocnseguently,  the  (ration)  supplies  are 
"pulled"  to  the  demand  locations  by  reguisitions.  During  this  process  the 
units  review  the  supply  reguirements  on  a  daily  basis,  if  possible; 
reguisitions  are  prs>ared  by  the  units  at  the  demand  locations;  the 
supplies  are  tran^xarted  (i.e.,  "pulled")  usually  by  truck  or  rail  finom 
the  storage  facilities  to  the  demand  locations. 

In  addition  to  the  push  and  pull  activities,  supplies  are  held  in 
inventory  to  ensure  preparedness,  preparedness  requiremmts 
determinations  are  based  on  troop  strmigths  f^*  a  given  level  and  duration 
of  conflict  and  current  feeding  policies  (such  as  MRE  to  T-Ration  ratios, 
etc.). 


Ration  Distribution  Structure 

The  structure  that  is  associated  with  rations  distribution  in  a 
theater  is  shown  in  Figure  1.  The  thin  arrows  represent  information  flow 
and  the  thick  arrows  represent  the  flow  of  (Class  1)  ration  supplies. 
Typically  the  information  flow  starts  at  tte  field  unit  and  reaches  TA  WC 
(Theater  Area  Materiel  Management  center) ,  which  is  the  integrated 
inventory  management  center  for  the  entire  theater.  O06OCM  WC  (Corps 
Support  OGnnand  Materiel  Management  Center)  is  the  functioned  control 
oenber  for  supply  and  malntenanoe  management;  it  relies  heavily  on  the 
ccBfxitBr  capabilities  of  GOGOCM  (Corps  Support  Ocnnand) .  C06OCM  ItC 
directs  storage  and  distribution  of  suteistenoe,  reviews  and  analyzes 
demands,  and  balances  the  workloads  of  subsistenoe  supply  units  at  or  near 
the  FICT.  ZMIC  (Division  Materiel  Management  center)  provides  centralized 
and  integrated  materiel  managaaent  support  for  Class  I  supplies.  OfC  is 
assigned  to  DUSOCH  (Division  Support  Conmand) .  The  WC  that  exists  in 
each  of  the  military  echelons,  i.e. ,  the  theater,  corps,  and  division 
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Theater 


Corps 


Division 


Brigade 


TA  MMC  =  Theater  Area  Materiel  Management  Center 
MMC  =  Materiel  Management  Center 
GS  =  General  Support 
DS  =  Direct  Support 

COSCOM  MMC  Corps  Support  Command  Materiel  Management 
Center 

DMMC  =  Division  Materiel  Management  Center 
DS  Main  =  Direct  Support  Main 
DS  Fwd  =  Direct  Support  Forward 


levels,  pcovides  its  oorres^nndijig  level  of  cxxjrdinatlon  of  simply 
activities.  Die  lower  (ndlitazy)  echelons  provide  Increasingly  more 
detailed  administration  of  supplies  in  their  own  areas.  The  IMC  of  each 
area  ocanunicates  with  Loth  General  Siqpport  (GS)  and  Direct  Sipport  (D6) 
units  and  provides  logistical  scpport  for  these  units.  Ohe  sipplies  axe 
physically  handled  only  by  the  GS  and  D6  units  (under  naraal 
circunstanoes) .  In  Figure  1,  note  that  D6  Main  provides  direct  scpport  to 
the  divisional  units  and  sends  siqpplies  to  D6  I\dd  (D6  Forward)  to  be 
distributed  to  all  units  at  and  below  the  brigade  echelon.  Also  note  that 
sipply  rcubes  that  are  exogenous  to  the  theater  ^pically  leexi  to  the  GS 
units  at  the  theater  and  corps  level;  these  supply  souxoes  can  originate 
fbOD  an  outside  nation  or  fxon  host  nation  sipport. 

The  generic  structure  of  a  distribution  node  is  shown  in  Figure  2.  As 
in  the  previous  figure,  information  flow  is  represented  by  a  thin  arrow 
and  (ra^on)  supply  flow  is  indicated  by  a  thiok  arrow.  Supplies  fron 
various  souroes  that  flow  into  the  area  are  inventoried  and  distributed  by 
the  GS  and  DS  units.  The  field  units  in  the  area  are  sipplied  by  the  DS 
units  and  report  supply  uszige  and  requisition  back  to  ttie  DS  unit.  Ihe  GS 
can  supply  the  DS  u^ts  in  the  area  and  the  GS  units  in  lower  echelons. 

Iha  generic  structure  of  each  location  is  represented  by  a  node  in  the 
logistics  network. 


Since  the  possibility  of  an  intense  nilltary  conflict  always  exists  at 
or  near  military  theaters,  the  ration  apply/distribution  system  must  be 
{dale  to  react  appropriately  to  provide  sudasistenoe  to  troops.  Iherefore, 
the  objective  of  the  logistics  nebwoiric  model  is  to  help  provide 
quantitative  decision  making  cipabllities  to  uppexHechelon  MC  users  in 
ttM  theater.  Ihe  top  priority  is  to  meet  the  demands  for  ration 
distribution  —  especially  the  demands  of  the  field  units  in  the  lower 
echelons 

Ihexe  are  two  major  ^pes  of  criteria  that  define  the  problematic 
issue  of  logistical  support:  prqaaredness  and  sustednability. 

Keparedness,  in  the  logistical  context,  is  the  ability  of  the  logistics 
system  to  react  effectively  and  decisiv^y  to  ary  adverse  cizcumstance 
that  may  occur  with  little  or  no  advance  warning.  Preparedness  includes 
the  ability  of  personnel  to  make  decisions  about  the  location  and  amount 
of  pre-positioni^  supplies,  and  the  ability  to  have  transportation  and 
disbribution  resources  in  place  and  in  order.  The  goal  of  preparedness  is 
to  provide  logistical  support  until  sustednability  cspebilities  axe 
achieved.  Sustainability  is  the  ability  of  the  logi^cs  tystem  to 
provide  long-term  supply  to  troops  once  the  %iartiffle  scenario  has  been 


aistaizmbility  {banning  sust  account  for  the  dynamic  and  probabilistic 
oexiditions  of  %i)ar.  The  unpredictable  characteristics  of  military 
conflicts  can  manifest  themselves  in  various  ways:  changes  in  intensity  or 
duration  of  battle,  (eneoy)  units  changing  loca^ons  rapidly  at  the  front 
lines,  consideration  of  rations  due  to  th^  vulnerability,  damage  to 
resouxoas  in  the  theater,  and  changes  in  demand  levels  due  to  troop 
casualties.  Characteristically,  the  lower  echelons  will  be  more  at  risk 
and  more  mobile  than  the  upper  echelons,  iihich  axe  relatively  stable  and 
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secure.  Sustednability  is  enhanced  vih^  the  ability  to  forecast  scenarios 
effectively  and  anticipate  troop  ncrvesoents  is  increased. 

Hie  human  element  can  undermine  sustainability  due  to  the  pursuit  of 
informal  activities,  such  as  trading  supplies,  or  sirply  by  ccraoitting 
(human)  errors. 

Another  factor  that  can  effect  ration  sustainability  adversely  is  the 
fact  that  wartijne  priorities  tend  to  favor  asnunition  and  fuel  over  ration 
supply  logistics. 

Model  Application  Strategy 

Tho^e  are  many  different  methods  that  are  useful  for  modelling 
aultiechelon  distribution  problems  that  have  tran^xirtation  lead  times  and 
inventory  levels  as  aspects  to  consider.  It  is  inportant  to  note  that 
althou^  there  are  many  kinds  of  anedytical  methods,  there  is  no 
established  method  that  can  cope  with  every  type  of  multiechelon 
inventary-dlstributlon  problem. 

One  analytical  technique  that  has  gedned  much  popularity  over  the  pem± 
dRcade  is  just-in-^ime  (JIT)  inventory  systems  analysis.  The  method  of 
JTT  is  inapplicable  to  the  ration  sipply/distribution  system  for  a  variety 
of  reasons.  First  of  all,  JTT  inventory  systems  are  strictly  "puli'* 
systems;  that  is,  they  respond  to  demand  and  they  assume  that  lead  time 
and  lead  time  variability  are  of  an  insignificant  magnitude  oonpared  to 
how  quickly  the  demands  can  change.  For  the  ration  si^ly/distslbution 
system,  that  assunption  does  not  epply  sinoe  lead  times  from  upper 
echelons  to  lower  echelons  can  take  days;  and  demands  at  the  front  lines 
can  change  much  more  reqpidly.  Secxsnd,  JTT  assumes  that  the  cost  of 
set*^,  ordering,  and  transpxjrting  are  insignificant  oonpared  to  the 
holding  oost  in  inventory;  again,  that  assuqjtion  clearly  does  not  apply 
to  the  ration  supply/distribution  system.  Even  though  JIT  analysis  has 
been  undergoing  some  mathematical  research,  it  is  not  yet  developed  encuc^ 
to  be  qpplied  to  the  logistics  problem  described  in  this  paper. 

The  modelling  technique  that  was  chosen  to  represent  the  theater 
logistics  is  a  minimum-cost  flow  networic  model  with  upper  and  lower  flow 
bound  representation  for  each  arc.  The  algorithm  that  is  used  as  part  of 
a  heuristic  to  reduce  the  operating  costs  is  known  as  the  out-of-ldlter 
algorithm.  Itoing  this  model,  inventory  distribution  facilities  are 
represented  by  nodes  and  transportation  routes  are  represented  by  arcs  in 
the  network.  To  take  the  time  element  of  the  logistics  system  into 
aooount,  the  method  of  dynamic  flow  was  applied  to  esqpand  the  model.  The 
details  of  the  logistics  network  model  are  described  in  the  next  section. 

An  overview  of  the  modelling  application  is  shown  in  Figure  3.  The 
schematic  illustrates  a  cyclical  process  in  which  information  about  the 
toeetor  (during  toe  military  conflict)  is  fed  to  the  database  in  the  form 
of  nahwcrk  parameter  vadues  that  define  the  current  situation.  The 
information  from  the  database  is  fed  to  toe  logistics  network  heuristic 
for  analysis.  Orders  and  instructions  that  are  based  on  toe  results  of 
toe  hsurlstic  anadysis  are  sent  to  toe  appropriate  parties  in  the  theater 
fkOB  toe  MNC  locations  in  the  echelons.  The  consequences  of  the  theater 
activities  are  aent  back  to  the  database  and  the  cycle  continues. 
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Figure  3. 

Modelling  Application  Schen\atic 
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lansncs  neiwsbk  mdcel 


Ite  sultiechelcn  logistics  model  is  depicted  mathanaticedly  as  a 
nsthPOKk  whose  basic  oomponents  consist  of  nodes  and  arcs.  Ohe  nodes 
retttssent  either  the  sipply  location  points,  such  as  Defense  logistics 
Acfency  (D[A)  depots  anVor  troop  issue  subsistence  activities  (TISAs) ,  or 
demand  points,  which  are  often  referred  to  as  ration  breakdown  locatim. 
The  arcs  are  vised  to  xepreBet±  the  sipply  routes  that  exist  between  the 
locations  (i.e.,  nodes)  in  the  theatre  of  operations.  Ihe  activity  that 
occurs  along  ttie  arc  that  connects  node  i  to  node  j,  that  is,  arc  (i,j)  is 
referred  to  as  the  flow  from  1  to  j .  Ihe  idea  is  to  perform  a 
■ultiohjective  optimlzatian  of  the  network  under  a  constrained  set  of 
■ultiooBDOdity  flows.  Each  ocmnodity,  k,  oocreqponds  to  a  sudo-network 
'Oiat  is  to  be  optimized  by  implication  of  a  minlirmn-oost-flow  algorithm. 


Ihe  primary  concern  of  the  model  is  to  ensure  that  the  proper  amount 
of  eaoh  ccnmodity  arrives  at  its  igpointed  destination  no  later  than  it  is 
demanded,  %henever  possible.  Ihe  svpply  and  demands  are  represented 
mathematically  by  constraints  within  the  network  model. 

Secondary  concerns  are  represented  by  the  multiok^ective  function  that 
exists  inherently  within  the  cost  stmcture  of  the  model.  Reducing 
various  kinds  of  product  flows  along  the  arcs  in  each  of  the  subnetworks 
enables  a  given  set  of  criteria  to  be  eoaphasized,  which  is  aooonplished  by 
manipulating  the  cost  structure  (of  the  arcs)  in  the  network. 

Ihe  four  criteria  of  fbndammttal  ixopcrtanoe  within  the  context  of  the 
avmply/distributlon  system  are  as  follows: 

(1)  Itanspoctation:  Reduction  of  the  strain  placed  on  the 
tranapartation  resources  (pertaining  to  truck  and  rail  oriented  srpply 
nechiffilsms)  is  iaportant.  Since  the  vehicular  resouroes  in  the  logisU.cs 
^^stsm  are  ^npically  the  most  critical  factor  in  the  operation  of  a 
theater,  this  crit^on  should  be  given  the  strongest,  or  nearly  the 
strongest,  priority  in  oonparison  to  the  other  criteria.  Placing  a 
relatively  high  set  of  (\itillty  function)  costs  on  the  arcs  that  represent 
transportation  flows  in  the  nebiork  will  place  enphasis  on  svppressing 
these  flows  in  oonparison  with  the  other  types  of  arc  flows  in  the 
network. 

(2)  Inventory  Levels:  Reduction  of  storage  in  pre  position  and  theater 
inventcndes  is  often  ueeful.  Reduction  in  overall  storage  levels  in 
theater  will  tend  to  speed  rotation  of  inventory  items.  An  incraase  in 
rotation  rates  is  partioularly  inportant  when  perishihle  items,  such  as 
rstions,  are  ixwolved. 

(3)  Demand:  Heeting  demands,  in  terms  of  preparedness  and 
sustainibility  can  be  aoocnplishad  by  effectively  distributing  the  stpply 
Ivmls  thioug^Kut  the  logistics  system.  Redistribution,  dLiring  each  time 
period,  in  reeponse  to  current  and  projected  demands  is  useful  in  avoiding 
stockouts,  lha  possibility  of  stockouts  is  a  particularly  critical  oonoern 
at  or  naar  the  front  line  of  troops  (FDGT) . 
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(4)  SiQjply  Production  Rates  and  lead  Tines:  Taking  sc^ly  souioe 
prodixition  rates  and  lead  times  into  oonsideration  is  very  important. 

Since  each  oonnodity  has  a  limited  rate  of  production  (in  OOMU5  plus 
OQCMUS) ,  demand  surges  can  overload  the  scqpply  souxoas  with  production 
demands.  If  the  ':pper  bounds  of  si:qpply  produi^on  rates  and  the  lead 

to  delivery  from  each  production  source  are  considered,  then  demands 
placed  directly  on  production  sources  can  be  smoothed  out  over  tine. 

Model  Sti\ifffarg 

Ihere  are  two  major  aspects  to  the  construct  of  the  logistics  network 
model.  First,  the  model  is  dynamic  over  a  countable  nuoher  of  tijoe 
periods;  therefore,  an  eaqpanded  version  of  ^le  network  representation  for 
one  time  period  is  used  to  model  the  system  over  the  entire  planning 
horizon.  Secsnd,  the  system  to  be  solved  is  deoonposed  into  \hat  will  be 
referred  to  as  ocnmodity  and  boundary  networks.  Ihe  volume-oonstrained 
boundary  network  takes  into  acoount  the  cunulative  effects  of  the 
volume-flow  fircm  all  the  ocnmodity  networks.  These  two  aspects,  referred 
to  above,  will  now  be  described  in  det2dl. 

Dynamic  Flows  Expansicn 

The  concept  of  the  dynamic  flow  version  of  the  (steady  state) 
logistics  network  model  is  more  easily  grasped  when  a  logistics  model 
having  only  one  ocnmodity  is  oonsider^.  Suppose  a  situation  exists  dose 
to  the  FLDT  in  which  a  small  section  of  the  logistics  ^stem  is  depicted 
by  Figure  4.  Nodes  nunhered  1  through  4  represent  storage  facilities. 

Node  S  r^sesents  an  exogenous  source;  that  is,  a  supply  source  that  mi^lht 
not  even  exist  in  the  theater  itself.  The  source  could  be  in  OMJS  or 
sonemhere  in  OCCMUS  near  the  theater.  Note  that  S  is  directly  supplying 
location  1,  in  Figure  4  as  shown  by  the  arc  that  connects  S  and  1.  Node  S 
does  not  have  to  represent  any  epedfic  location;  the  idea  is  that 
location  1,  in  Figure  4,  is  being  supplied  by  a  source  or  gro^  of  source 
locations  that  are  represented  by  node  S;  arc  (S,l)  conveys  ttiis  concept. 
Nodes  P  and  D  denote  pre  position  and  demand  at^'vdties  respectively. 
Neither  P  nor  D  represent  any  location  at  all;  instead,  the  arcs  that 
emanate  from  node  P  and  whatever  arcs  that  enter  node  D  represent  supply 
and  demand  flow  activities,  respectively.  The  activities  of  all  thB  arcs 
axe  constrained  by  upper  and  lower  bounds  and  regulated  by  costs;  this 
will  be  described  later  in  this  chapter.  The  lunbers  on  each  arc  of 
Figure  4  denote  the  nanher  of  time  periods  xeguired  for  a  supply  shipment 
to  travel  along  the  arc  tvaa  node  to  node.  Seme  of  the  arcs  axe  assigned 
(time)  values  of  zero  because  the  activity  that  the  arc  represents  tal^ 
no  significant  time  to  occur.  For  instance,  arc  (4,D)  xegoixes  no 
significant  time  since  it  does  not  represent  an  actual  route  in  the 
logistlcB  eystem.  Arc  (4,D)  exists  as  an  activity  to  ensure  that  the 
total  demand  of  the  oonmodity  flawing  into  node  (location)  4  from  nodes  2 
and  3  does  not  fall  beneath  a  prescribed  lower  limit  (that  is  given  by  the 
lower  bound  on  arc  (4,D) ) . 


II 
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Figure  4. 

Network  Representation  of  Four  Facilities  Near  the 
Front  Line  of  Troops  with  their  Supply  Routes 

and  Activities 

Nodes  1  through  4  represent  storage  facilities. 


Ihe  dynamic  flew  verslc3n  of  the  graph  in  Figure  4  is  shown  in  Figure 
5.  Note  ^lat  there  are  nodes  that  han«  two  nunhers  in  their  designations; 
the  first  nuoiaer  denotes  the  tine  period  and  the  seoend  nunber  corresponds 
to  Idle  Identifying  nunber  of  the  node  in  the  grs^  in  Figure  4.  Ihe  P,  S, 
and  D  nodes  denote  the  same  meanings  as  in  Figure  4.  The  node  referred  to 
in  Figure  5  as  SS  is  the  si^ersouroe  node;  although  SS  does  not  represent 
any  real  location  the  arcs  emanating  from  it  have  a  fUnctionad.  meening. 

Ihe  only  reason  node  SS  is  important  is  due  to  the  fact  that  network 
algorithms  typically  solve  graphical  representation  of  networks  having  one 
souroe  node,  namely  SS,  and  one  terminal  node,  namely  D.  Networks  having 
aultiple  sources  and  multiple  terminals  (or  siidcs)  are  typically  modelled 
using  siperscuroe  and  sipersihk  nodes  in  the  2bove  manner. 

Ihe  oonoq±  of  the  eopanded  graph  (see  Figure  5)  can  be  grasped  by 
considering  seme  brief  exanples.  For  instance,  consider  nodes  l  and  2  in 
Figure  4.  Since  one  time  period  is  required  to  send  a  shipoaent  over  arc 
(1,2) ,  this  actively  is  rqmesented  in  Figure  5  by  arcs  fixm  node  (1,1)  to 
node  (2,2),  from  (2,1)  to  (3,2),...,  ton  node  (t,l)  to  node  (t+1,2), 
etc.,  from  node  (5,1)  to  node  (6,2).  In  another  example,  consider  nodes  2 
and  3  in  Figure  4.  Sinoe  two  time  periods  axe  required  to  send  the 
oemmodity  along  arc  (2,3) ,  the  r^mW^ntation  on  the  esqpanded  grsqph  are 
arc  from  nodes  (t,2)  to  nodes  (tf2,3)  fear  t  «  1,  2,  3,  and  4  in  Figure  5. 
In  general,  if  arc  (A,B)  from  node  A  to  node  B  requires  h  time  periods  for 
a  oenmodi^  to  traverse  it,  then  the  representation  %Kuld  be  an  arc  ton 
node  (t,A)  to  node  (bth,B) .  Over  a  planning  horizon  of  T|{,  that  is, 
time  peil.ods  1,2,...,  T|],  t  would  be  defined  as  an  integer  such  that 
1  ^  t  <  Tji-h.  In  Figure  5,  (for  i  -  1,...,  6)  is  the  demand  for 
the  oonmndlty  at  time  period  i. 

Ihe  logistics  network  model  takes  into  consideration  five  types  of 
activities;  th^  are  described  as  follows  (using  the  notation  of  the 
dynamic  flows  greph) : 

(1)  Transportation:  Arc  ((t,A),  (t4h,B))  as  described  ebove. 

(2)  Storage:  Arc  ( (t,A) ,  (t'fl,A) )  for  each  time  period  for 
t  *  1, . . . ,  %-l,  for  location  (node)  A. 

(3)  Demand:  Arc  ( (t,A) ,  D)  for  demand  at  node  A  for  t  »  1, . . . ,  Tjj. 

(4)  Sources:  Arc  (S,  (t4h,A) )  for  ahipnents  ton  the  exogenous  source 
to  node  A  for  t  >*  1, . . . ,  Tjj-h. 

(5)  lie-positioning:  Arc  [P,  (1,A)  ]  for  all  locations,  node  A,  where 
pre  positioning  is  oensidered. 
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SS  =  Supersource 
S  =  Exogenous  Source 
P  =  Pre-position 
D  =  Demand 


Figure  5. 

Dynamic  Flows  Expansion  of  the  Network 
in  Figure  4  for  Six  Time  Periods 

In  the  intermediate  nodes,  the  first  indice 
represents  the  time  period  and  the  second 
number  represents  the  identifying  number 
of  the  node  from  Figure  4. 


So  far,  the  axes  involved  have  been  described  as  routes  that  require  h 
tine  periods  to  traverse.  There  are  additional  paranebers  that  define  the 
arcs  that  exist  nathenatically  in  the  model.  The  parameters  of  arc  (i,j) 
for  the  (eaqpanded)  dynamic  flows  graph  are  as  follows,  for  each  oenmodity 
k: 


(1)  The  lower  bound  of  the  flow  is  l(i, j,k) . 

(2)  The  upper  bound  of  the  flow  is  u(i,j,k). 

(3)  The  cost  per  unit  of  the  flow  is  c(i,j,k}. 

The  quantity  f(i,j,k)  is  the  nunher  of  items  (e.g.,  xatiens)  of 
oonnodity  k  thet  flow  along  arc  (i,  j)  during  the  time  periods  that  are 
represented  by  arc  (i,  j) .  It  follcws  that  the  flow  must  satisty  the 
condition: 


<  f(i,j,k)  <  u(i,j,k)  (1) 

for  2dl  arcs  (i,j)  for  every  ooenodily  k.  The  upper  bound  (e.g., 
u(i,j,k))  is  usually  referred  to  as  "ceqpacity"  of  the  arc. 

The  cost  of  the  flow  of  arc  (i,j,k)  is  equal  to  c(i,j,k)  f(i,j,k). 
The  objective  is  to  minimize  the  sum  of  edl  the  arcs  for  each  oonm^ty 
k  *  1, . . . , 


Minimize  S  £  c(i,j,k)  f(i,j,k)  , 

i  j 


for  all  k.  (2) 


The  parameters  that  are  assigned  to  each  arc,  in  the  dynamic  flows 
model,  distinguish  the  activity  of  each  of  the  arcs;  the  manner  in  which 
this  is  done  is  described  in  the  following  five  activities: 

(1)  For  tran^nrtation  activities,  the  lower  and  upper  bounds  are  zero 
and  infinity  respectively.  The  actual  nunber  of  each  oenmodity,  k,  being 
transported,  in  itself,  is  not  a  limiting  factor;  the  volume  (or  cube)  is 
a  pcoblsnatic  aspect  in  the  logistics  system.  There  are  wit  costs  that 
can  be  assigned  to  each  oenmodity  to  discourage  the  use  of  tran^artatloo 
rescurns.  The  unit  cost  c(l,j,k)  can  be  bas^  on  the  mass  and  volume  of 
oenmodity  k  as  well  as  on  considerations  that  relate  to  the  particulars  of 
the  route  of  arc  (i,j)  and  any  utility  functions  that  the  users  of  the 
model  dean  fit  to  construct. 


(2)  For  storage  activities,  the  lower  bound  is  set  to  zero,  and  the 
upper  bound  u(i,j,k)  is  typically  set  to  infinity,  wiless  other 
considerations  that  axe  inherent  to  the  paxtic^ar  kind  of  oenmodity  are 
relevant.  As  in  the  tionsportation  activity,  it  is  the  unit  volume  of 
each  oenmodity  that  relates  to  the  nunber  of  items  that  can  be  stored  at 
location  i  from  one  time  period  t  to  the  next  tfl.  Unit  costs  can  be 
assigned  to  ttie  arcs  to  influence  the  amount  of  items  being  stared.  These 
costs  can  be  related  to  a  utility  function  that  is  selected  by  the  users 
of  the  model. 
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(3)  For  denentl  activities,  the  arcs  have  a  finite  positive  lower 
bound,  recall  dj^  in  Figure  5  for  the  beginning  of  time  period  i. 

Setting  the  parameter  l^i,j,k)  equal  to  the  demand  for  ocninodity  k  at  node 
i,  tdieire  j  corresponds  to  node  D  in  Figure  5,  forces  that  demand  to  be 
met.  ISie  c^eKdty  u(i,j,k)  is  infinite.  The  cost  is  set  ec^ial  to  zero. 
The  demand  is  tdiat  drives  the  flow  thrcu^  the  system.  Note  that  in  the 
mazketplaoe,  where  these  ocnnodities  are  sold,  negative  costs  would  be 
assign^  to  enoourage  sales. 

(4)  For  activities  that  r^jresert  sources,  a  lower  bound  of  zero  and  a 
finite  positive  vpper  bound  is  assigned,  that  is,  u(i,j,k)  >  0.  The  iqpper 
bound  is  the  capacity  of  the  souroe  to  produoe  flov,  f (i,j,k) ,  vdiere  i 
oorreeponds  to  the  source  node,  S,  j  oorreqponds  to  the  receiving  facility 
(node) ,  and  k  oorreaponds  to  the  oonnodity.  The  unit  cost,  c(i,  j,k) ,  can 
be  set  to  relate  to  the  market  prices.  Note  that  in  a  maj(»:  wartime 
eaergency,  the  costs  may  be  assigned  zero,  since  the  lives  of  the  soldiers 
and  the  security  of  regions  that  are  affected  most  direc±ly  oculd  make 
capitalistic  issues  insignificant  in  oocparison. 

(5)  For  pre-positioning  activities,  the  arcs  would  be  assigned  bound 
parameter  values  that  are  analcsgous  to  the  values  of  the  source  activity 
arcs,  that  is,  a  lower  bound  of  zero  and  an  vpper  bound  u(i,j,k)  >  0. 

Node  i  oorzes|X3nds  to  node  P  and  node  j  oorr^ponds  to  the  facility  having 
pone-positioned  u(i,j,k)  units  of  ocnnodity  k  just  prior  to  the  beginning 
of  the  first  time  p^od.  No  casts  would  be  assigned  the  activity;  that 
is,  c(i,j,k)  »  0,  unless  the  effcn±  to  remove  the  items  from  a 

pre-  positioning  mode  of  storage  to  an  orientation  that  cxarresponds  to 
storage  that  has  a  higher  level  of  readiness  is  considered  costly  for  seme 
particular  reason. 

It  is  inportant  to  note  that  arcs  (SS,P)  and  (SS,S)  in  Figure  5  have 
lower  and  rpper  bounds  of  zero  and  infinity,  respectively.  Also  note  that 
the  unit  cast  of  both  the  arcs  is  zero. 

Note  that  there  could  be  five  types  of  exsts  that  are  asscaiated  with 
the  loglstias  network  model;  each  type  of  aoet  corre^ands  to  its 
asscaiated  activity.  Because  of  the  different  types  of  casts  in  the 
network  model,  this  model  can  be  classified  as  a  multlcli^ective 
optimization  model  in  which  each  objective  oorre^ands  to  one  of  the  five 
types  of  activity  casts.  As  e>g)ledned  in  the  introduction,  the 
predoninant  costs  in  the  (military)  theater  would  be  transportation  and 
(inventory)  storage  costs,  na-posltioning  costs  may  be  a  significant 
factor  in  acme  theater  scanarics. 

Effects  af  Qgimfrfttv  Vbluroe-Flcx» 


A  major  aspect  to  the  structure  of  the  model  is  the  cumulative 
influence  that  the  vcdune-flaw  of  all  ocemodities  has  with  respect  to  the 
vcdUMtrlc  capacities  of  each  of  the  routes  in  the  logistics  network.  In 
eesenoe,  the  issue  is  cubage.  Each  oconedity  has  a  oorre^xnding  volume 
(or  odbsge  V-(k) ,  for  ooewodities  k  •  1, . . . ,  N^,  where  N^  is  the 
number  of  different  cxnmodities  in  the  logistics  tystan.  The  cumulative 
VQlune-flcw  through  each  arc  (i,  j)  in  the  tynamic  flows  network,  when 
sumwed  over  all  the  ooenoditim,  is  given  by  equation  (3) ;  such  a  network 
of  omulative  volume-flows  is  referred  to  as  a  boundary  (flow)  network. 
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(3) 


v(i,j) 


c 

2  f(i,j,k)  Vc(k) 
te=l 


Recall  that  each  oomoodity,  k,  has  its  cwn  oarxe^xnUjig  dynamic  flows 
network.  Mhen  each  of  the  flows  f(i,j,k)  in  eOl  dynamic  flows 
networks  is  determined,  the  cunulative  vol\ane>flows  for  all  arcs  (l,j)  in 
the  dynamic  flows  (boundary)  network  axe  ocofubed  and  ocnpaxed  to  their 
oorreqpcnding  volumetric  flew  capacities,  U^(i,  j) ,  in  the  dynamic  flows 
(boundiuy)  network  for  (cunulative)  volume-uows. 

It  is  iirportant  to  note  that  the  dynamic  flows  network  for  volumetric 
flows  has  exactly  the  same  graphical  structure  as  each  of  the  graphs  for 
the  dynamic  flows  of  the  oaraaodities. 

It  is  necessary  to  ensure  that  none  of  the  voliaaetric  flows,  v(l,j), 
exooods  the  vpper  bound,  Uy(i,j),  in  the  volumetric  dynamic  flows 
network.  Ihe  heuristic  for  nulticcnnodity  flows,  ihich  is  described  in  a 
later  section,  identifies  ary  violation,  such  that  v(i,  j)  >  U^(i, j) ,  and 
eliminates  it  during  the  next  Iteration  in  the  heuristic. 

nixee  of  the  five  types  of  activities  in  the  logistics  network  have 
finite  volumetric  flow  capacities;  they  are  transportation,  storage,  and 
pre-positioning.  Since  demand  and  source  activities  do  not  have  physical 
routes  tiiat  axe  constrained  by  ary  practical  considerations  with  reepect 
to  volume  (i.e. ,  cubage) ,  the  arcs  that  oarxespond  to  these  activities  are 
oonsidered  to  be  effectively  boundless.  Thus,  U^(i,  j)  is  set  to  infinity 
(or  an  extremely  large  positive  quantity  that  could  not  possibly  be 
reached  by  the  volume-flows) . 


3iQ2lu 


Flew  TAmii-Jit-AnnR 


The  network  flow  activities  that  have  finite  volumetric  cspaclties 
have  their  oarte^)onding  upper  bounds  for  various  reasons  and  with 
differing  characteristics.  For  instance,  storage  and  pre-positioning 
limits  are  merely  static  in  oenparison  to  tran^xsrtation  since  the  nunter 
of  vdiicles  at  a  facility  (location)  can  change  very  rapidly  fixm  cne  dty 
to  the  next.  Storage  and  pre-positioning  (volumetric)  cepacities  have 
essentially  the  same  characteristics  as  each  other;  they  most  strongly 
xelatte  to  the  internal  dimensiens  (i.e.,  actual  size)  of  the  storage 
fedlity  that  determines  their  V2dues.  For  that  reasen,  the  volumetric 
capacities  associated  with  storage  ferdlities  are  nearly  constant,  but  not 
always  constant  sinoe  catastrophic  events  on  the  battlefield  can  suddenly 
reduce  ttieae  capacities  to  zero  or  near  zero.  If  eapected  capacities  are 
weed  when  it  is  decided  that  the  ri^  (or  probability)  of  destruction  is 
significant,  then  these  values  will  indeed  tend  to  fluctuate  more  actively 
from  one  time  period  to  the  next. 
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Transportation  (volumetric)  cspacities  are  a  much  more  elaborate  issue 
than  those  of  storage  and  pre-positioning,  consider  a  route  frcn  location 
i  to  destination  j,  within  the  context  shown  in  Figure  4.  Ihe 
periods  required  for  the  vehicles  to  get  frxxa  i  to  j  and  return  to  i  are 

denoted  and  ,  re^ectively.  The  siperscripts  D  and  R  stand 

for  destination  and  return,  reflectively.  T|^j  is  the  quantity  that 

affects  the  time  dynamic  in  the  dynamic  flows  network.  Note  that  h,  «hich 
was  defined  in  a  previous  section  as  the  number  of  periods  required 

to  send  a  oomnodity  along  arc  (i,  j) ,  has  the  same  meaning  their 

oontexts  obviously  differ  from  one  another.  Since  the  v^des  (usually 
tnicJcs)  are  not,  as  a  rule,  lent  out  to  other  depots,  they  return 

(T^^j  +  T^^j)  time  periods  erfter  leaving  location  i.  Denote 

fR) 

as  the  transportable  voltmie  cepacity  at  location  i,  and  vl  4  as  the 

J 

volume  that  is  actually  transported  &cin  i  to  j;  ecudi  starting  at  the 

fR) 

beginning  of  time  period  t.  Note  that  is  not  less  than  the 

actual  volume  of  products  being  ^pped  from  facility  (location)  i  to  j 
during  t.  Ihus,  location  i  will  lose  ^  volumetric 

(transport)  capacity  at  time  period  tfl  would  be: 


„(C)  .  „(C)  (R) 

Vi,t+X  - 


M) 


if  the  volunetric  gain  from  the  returning  vehicles  from  location  j  were 
not  considered.  Ihe  volximetric  gain  fron  the  returning  vdiides  is 

(R) 

denoted  by  j  Subscript  q(i,j)  is  defined  as 

q(i, j)  -  (t+1)  -  (T^^J  +  T^^^)  .  (5) 

Note  that  q(i,  j)  must  be  a  positive  integer.  Iherefore,  taking  the 
rsbuming  vehicles  into  account,  the  volumetric  cepacity  is; 


^(C)  ,  „(C)  _  (R)  (R) 

Vi,t+1  -  %t  -  Vi,j,t  +  ^i,j,q(i,j)* 


(6) 
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Iflcation  1  at  tine  t  and  sixxsessive  time  periods,  tfl, . . . ,  Tu  ate 
related  to  the  stxvictuie  of  the  dynamic  flows  network  as  described  in  a 
previous  section  (of  this  pe^)er)  that  discusses  the  structure  of  the 

(C) 

networks.  Therefore,  the  calculated  values  of  for  all  i  and  t, 

can  be  assigned  to  the  corresponding  czpacities  U»(i,j)  of  the  (^namic 
flows  network. 

As  previously  noted,  tAienever  the  volisnetric  flow,  v(l,  j) ,  of  the 
dynamic  flows  netwoidc,  exceeds  U^(i,j),  corrective  measures  are 
i^lenented  iteratively  by  using  the  heuristic  for  iwii'UnrBMnrri'i't-y  flows. 
The  heuristic  is  described  in  detail  in  the  next  main  section. 

Solution  Scheme 

canoe  the  ooranodity  flow  configuration,  f(i,j,k),  for  all  arcs  (i,j) 
and  oonnodities  k  *  1, . . . ,  in  the  dynamic  flows  network  mnAai  are 
determined,  reoonnendatlons  are  made  in  terms  of  orders  and  instructioi^; 
refer  to  Figure  3.  The  reooDinendations  are  based  on  the  solution  and  can 
be  qualitative  as  well  as  quantitative.  The  reconmandatiots  are  received 
by  theater  logistics  personnel  and  are  isplenented,  either  fully  or 
partially,  according  to  the  judgement  of  these  personnel. 

Wien  the  next  time  period  arrives,  the  current  status  of  the  logistics 
netenrk  is  fed  back  to  the  logistics  database  and  the  current  status  of 
the  distribution  supply  eystan  is  used  to  i:pdate  the  parameter  V2dues  of 
=he  dynamic  flews  network  model.  Agzdn,  the  flow  configurations, 
f (i, j,k) ,  are  oonputed;  and  the  cycle  continues  as  long  as  needed,  as 
depicted  in  Figure  3. 


i 
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MJLnCXlfODm  NEIWORK  HEURISTIC 


The  heuristic,  vftilch  is  outlined  in  this  repeat,  uses  a  network 
optimization  technique  that  is  known  as  the  out-of-kilter  edgorithm^'^ 
as  part  of  its  solution  methodology.  Ihe  algorithm,  which  will  be 
referred  to  as  OKA  for  brevity,  is  useful  in  optimizing  min^mim  cost  flow 
network  models.  Each  oonmodity  in  a  oultiocoiiKxiity  network  is  solved 
separately  in  the  heuristic.  Ihe  rest  of  the  methodology  in  the  heuristic 
involves  the  elimination  of  constraint  violations  «hen  the  volumetric  flow 
v(i,j)  exoeodfl  the  volumetric  flow  cepacity  U^(l,j)  of  arc  (l,j).  If 
tl^  violations  cannot  be  removed  zifter  a  predetermined  nunber  of  attaipts 
(i.e. ,  cycles)  or  the  oonmodity  network  themselves  are  infeasible 
initiidly,  the  lieuristic  will  terminate  without  a  final  solution  to  the 
ooHiBfxIity  network  flows,  f (i, j,k) . 


The  out-of-kilter  algorithm  (0K2V)  is  initiated  with  a  starting 
feasible  solution  in  the  (primal)  flow  variables.  The  dual  variables  are 
then  calculated  and  the  constraint  violations  in  the  dual  constraint  set 
are  assess^.  The  assessment  is  made  using  quantities  called  kilter 
nuBtars.^'^  Each  arc  (i,  j)  in  the  graph  has  its  cDtrespcnding  kilter 
nuDber.  Mien  the  flow  in  the  network  is  augmwited  81.J1  that  ^1  the 
kilter  nuDbers  are  zero,  dual  feasibility  is  achieved;  and,  thus,  the 
optimal  (flow)  solution  is  determined.  A  set  of  dual  vari^les,  Miose 
values  are  su^  that  each  variable  oorrespends  to  each  node,  uniquely 
determines  the  optimizing  solution  (i.e. ,  analogous  to  fingerprinting  or 
ENA  coding  in  determining  each  individual's  identity) .  In  performing  the 
iterations  with  GEKA  it  is  found,  using  the  (conputer)  code  that  was 
written  as  part  of  the  project,  that  Mien  the  largest  kilter  nunber  is 
selected  to  be  reduced  to  zero  the  algorithm  oenverges  significantly  more 
quickly  than  Mien  the  minimum  kilter  nunber  is  chosen  to  be  reduced  to 
zero.  Note  that  kilter  nunbers  can  never  have  negative  values. 

The  advantages  of  using  OKA  are  considerable.  Pertiaps  the  most 
inportant  reason  to  oonsider  the  use  of  OKA  is  that  lower  (flow)  Jxunds 
can  be  defined  on  each  arc;  with  the  mininum-oost-flow  algorithiD^  the 
modeller  has  to  settle  for  zero  as  the  lower  bound.  Another  useful 
characteristic  of  OKA  is  that  noninteger  costs  (i.e.,  c(i,j)  fren  the 
previous  chapter)  can  be  wed  without  jeopardizing  oonvergeixx  to  the 
optiaun  after  a  finite  nmber  of  iterations.  Another  advantaige  in  wing 
OKA  is  that  when  a  change  in  a  parameter  vadw  (such  as  a  bound  or  a  unit 
cost)  has  been  made  after  the  solution  of  the  original  network  has  been 
obtadned,  the  dual  variable  values  that  oorrespond  to  the  nodes  in  the 
network  oeoi  be  used  as  a  starting  solution  for  the  altered  network.  Ey 
applying  the  dual  variable  values  to  the  altered  network,  a 
charmctaristioally  rapid  oonveegenoe  is  achieved  to  obtadn  the  new  optimum 
doer  oenfigparation,  that  is,  previdsd  the  changes  did  not  render  the 
altered  nebwerk  InfeasiliLe. 
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HaiiHgfcle  Sta-niriteune 


The  a$)plication  of  to  the  solution  of  a  oultiocninodity  network 
flow  prcblen  is  outlined  in  detedl  in  this  section.  First,  it  should  be 
noted  that  although  there  are  other  solution  techniques  that  solve  scoe 
types  of  miltioonnodlty  flow  problans,  they  are  limited  lainly  to  the  sums 
of  the  flows  themselves,  that  is,  f(i,j,l)  f(i,j,2)  f(i,j,NQ) 

and  that  the  methodologies  tend  to  be  centered  on  sii^ect  areas  covered  by 
linaar,  integer,  and  large-scale  progranning  methods.  The  following 
heuris^  entails  a  much  different  i^pproach  than  the  mathematical 
progranning  methods  and  converges  to  a  solution  relatively  quickly  sinoe 
it  usee  OKA  as  a  basis  for  its  seardh  routine. 

The  heuristic  will  be  described  here;  Figure  6  can  be  used  as  a 
guide.  The  main  idea  of  the  heuristic  is  to  invoke  a  new  and  possibly 
tenporary  constraint  to  the  arc  that  esdhibits  the  maximum  violation  of  its 
v^pper  (voliimetric)  cE^iacity  bound;  recall  11^(1,  j) .  Note  that  maximum 
(constraint)  violations  were  tested  in  the  heuristic  with  a  rate  of 
convergence  that  was  significantly  faster  than  in  the  case  of  correcting 
violations.  The  added  constraint  %i»s  in  the  form  of  a  change  in 
one  or  more  of  the  vpper  bounds  in  the  oonnodity  networks;  it  must  be 
that  the  added  constraints  were  labeled  as  not  being  part  of 
tile  original  formulation  with  respect  to  tiie  bound  parameter  values. 

These  a^usted  bounds  were  determined  by  using  a  projection  method;  the 
projection  method  is  described  in  detail  in  the  i^jpendix.  After  the 
adjusted  UEfier  bounds  are  added  to  the  oomnodity  networks,  k,  for 
k  »  1, . . . ,  Ng,  OKA  is  iaplemented  agsdn  and  the  cycle  continues,  if 
there  are  no  other  oonplicatlons.  There  can  be  many  ocmplications, 
however,  in  a  typical  problem  of  this  kind. 

The  most  basic  of  the  ocoplications  referred  to  above  is  the 
occurrence  of  too  many  bound  adjustments  added  to  the  commodity  networks 
resulting  in  the  flows  being  choked  off  from  other  arcs  in  tiie  networks. 
Older  these  conditions,  the  demands  that  are  inherent  from  the  lower 
(flow)  bounds  of  some  of  the  arcs  cannot  be  met;  therefore,  infeasibility 
results  in  one  or  more  of  the  over-constrained  ocnmodity  netwoadcs.  To 
remedy  that  situation,  all  ipper  bound  (constraint)  adjustments  prior  to 
the  most  recent  adjustment  are  deleted.  The  bounds  on  the  remaining 
adjusted  arc  are  then  "fixed";  this  means  that  the  lower  bound  (in  each 
ocnmodity)  is  set  equal  to  the  adjimted  ipper  bound.  This  will  fix  the 
amount  of  flow  through  the  arc  to  a  oenstant  value.  Note  that  the  term 
"constraint"  in  Figure  6  and  the  terms  "bound  CK^ustment"  axe 
interchangeedale  in  the  description  of  the  heuristic. 

The  deletion  of  all  but  the  most  recent  bound  ai^ustment  mey  hzqppen 
more  than  onoe  during  the  heuristic's  search;  therefore,  the  erijusted  or 
"fixad"  arc  that  is  left  undeleted  eifter  all  the  other  bound  adjustments 
axe  erased  mey  be  the  same  (remaining)  adjusted  arc  as  earlier.  This  is 
"cycling."  The  identity  of  all  the  fixed  arcs  is  recorded  to 
detect  eyeing. 
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Figure  6, 

Multicommodity  Network  Heuristic  Overview 
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Itien  an  arc  is  identified  as  a  "cycling  arc”,  its  bounds  (upper  and 
lower)  ronain  fixed  at  the  sane  value  and  the  arc  is  "partitioned"  away 
fron  its  original  set  of  arcs.  Ohe  partitioning  neans  that  the  arc  joins 
a  collection  (i.e. ,  another  set)  of  arcs  that  will  never  be  altered  in  any 
way  as  long  as  the  heuristic  is  en^tged  in  solving  that  particular  network 
problen.  If  the  nuiDber  of  partitioned  arcs  exceeds  a  prescribed  limit, 
the  heuristic's  search  is  terminated.  A  raixmiiended  value  for  the  upper 
bound  for  cycling  arcs  is  three,  which  is  based  solely  on  test  nsis  that 
wiere  perfoKBed.  Several  e}g)erimental  network  prcblens  that  were  very 
tightly  constrained  yet  known  to  be  feasible  were  perfomed;  all  of  these 
networks  yielded  feasible  solutions  in  two  or  less  discoveries  of  a 
cycling  arc. 

Solution  Characteristics 

Sixne  each  of  the  ccratcdlty  networks  is  optimized  separately  and  under 
the  presence  of  vpper  and  lower  bound  adjustments,  the  solutions  to  the 
BUltioomnodity  networks  must  be  considered  subcptimal  in  general. 
I%eliminary  trials  with  some  network  problems  yield  results  that  suggest 
that  the  heuristic  frequently  comes  i^thin  one  percent  of  the  true  optima 
even  though  many  of  the  esqierimental  problems  were  very  tightly 
constrained. 

Data  Manipulation  for  Meraorv  Oonsiderations 

Since  one  of  the  goals  of  the  project  was  to  allow  the  software  to  be 
used  in  many  different  circumstances  and  at  many  echelons,  an  attenpt  was 
made  to  make  the  software  practical  for  use  on  a  personal  oonputer;  for 
this  reason,  data  files  were  manipulated  by  the  software  in  order  to  save 
apace  in  storage  within  the  executable  code  itself. 

In  qpite  of  the  above  ektjustments,  some  theater  logistics  models  may 
have  too  many  facilities  for  the  software  to  deal  with  cn  a  p**'"'*'*^'* 
oonputer.  It  may  be  useful  to  solve  the  dynamic  flows  network  model  in 
sec^ons,  under  the  conditions  just  mentioned.  Ihe  solution  at  the 
boundary  to  a  given  section  mey  be  used  as  an  input  to  the  boundary  of  an 
adjacent  section.  lunping  facilities  that  are  located  relatively  close  to 
each  other  in  the  theater  may  be  another  effective  method  in  circumventing 
such  a  problematic  issue  when  a  personal  oonputer  is  eaployed. 
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CCNODSICKS  AND  REQdfOlDATIONS 


Althou^  a  generic  model  in  the  form  of  a  network  vas  defined 
oatheBBtic^y  and  a  heuristic  was  developed  that  is  ce^mble  of  producing 
a  adixtion  to  the  laiLtiocninodity  flow  nebKadc,  further  esqperimenibtticn  of 
the  model  and  heuristic  va^  yield  additional  matheDatical  insights, 
currently,  further  refinement  of  the  model  and  heuristic  is  not  a  critical 
issue.  The  {plication  and  extension  of  the  model  and  heuristic  software 
to  either  an  actual  or  a  small  trial  problem,  «hich  is  similar  to  an 
actual  or  a  historical  theater  logistic  problem,  are  the  next  inportant 
steps.  Further  development  would  lead  toward  an  application  to  an 
existing  thmiter. 

Further  reoonnendations  can  be  made  to  enhance  the  effectiveness  of 
the  model  fay  adding  technical  extensions.  For  instance,  the  development 
or  ap|d.ica1^on  of  forecasting  technigues  that  use  artificial  intelligence 
as  a  major  part  of  their  methodology  may  provide  interesting  and  useful 
insights  and  ce^iabilities.  In  a  similar  vein,  a  set  of  "learning” 
parameters  that  relate  to  the  suooeas/fedlxxre  history  of  the  logistics 
model  could  possibly  be  developed.  The  values  of  the  learning  parameters 
could  be  altered  heuristicEdly  to  influence  the  future  response  of  the 
model. 

The  capabili^  of  the  model  could  be  extended  if  input  from  domain 
experts  is  xised.  Domain  experts  are  individuals  iho  possess  a  deep 
tmderstanding  of  the  oonplexltles  and  intricacies  of  a  particular  subject 
area.  Theater  IfC  personnel  and  inventory  managers  may  be  enployed  as 
domain  eoperts  for  the  model;  their  input  may  be  added  to  the  model  in  the 
form  of  syllogistic  rules  and  guidelini^  to  extend  the  existing  rules  base. 

If  the  model  and  its  related  software  are  enhanced  by  sene  or  all  of 
the  ideas  that  are  described  above,  then  it  oould  also  have  high  civilian 
potential  as  a  useful  tool  in  deeding  with  inventory  and  transih^nent 
problems  in  the  private  sector. 
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APPENDIX 

PROJECnON  METIHOD 
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APPENDIX 


FPCOECnCH  MEIKX) 


Ihe  projec±ion  xnethod  is  used  to  ai^ust  the  upper  bounds  of  arcs  in 
th&  oQwnodity  networks  so  that  the  volumetric  flow  v(i,  j)  is  lowered  to 
e(]oal  the  volumetric  flow  ce^iacity  Uy(i,j)  of  arc  (i,j)  in  the  dynamic 
flows  network. 


For  oonveniencse  in  notation,  denote  vector  fj^^ j  as 

f (i» j »2) , . . . ,  f(i,^,NQ))  ,  (A“l) 

where  is  the  nmober  of  ocmnodities  in  the  netwradc  system.  Similarly 
denote  as 

(V^jCl),...,  Vj,(Nc))  .  (A-2) 


Also,  denote  U^(i,j)  as  Uj^^j. 

Hie  idea  is  to  project  vector  f  j  in  the  direction  of  negative 
fay  a  euclidean  distance,  r,  such  that*' 


(A-3) 


Hiis  is  referred  to  as  (the)  projection,  and  f^^  j  is  referred  to  as  the 

adjusted  vector.  Note  that  the  dot  product  of  f|^i  and  has  the 
same  meaning  as  equation  (3)  in  the  section  that  describes  the  logistics 
network  model.  A  conplicaticn  arises  because  none  of  the  flow  oonponents 
is  allowed  to  become  negative.  Because  of  the  non-negativity  resection, 
the  heuristic  farces  some  of  the  conpcnents  (i.e. ,  the  ccnaodities)  of  the 
vector  to  remain  static  ^Aiile  other  oonponents  are  allowed  to  vary  as 


part  of  the  projection. 


Define  f|^j 


as  shown  below: 


— (v) 

f^^j  *  (f  {i»  j » 1) » •  •  • »  0,...,  0)  (A— 4) 

and 

"(S) 

H,j  -  (0,...,  0,  f(i,j,iiri-l),...,  f{i,j,Nc)),  (A-5) 


in  which  the  first  m  ocaponents  are  permitted  to  vary  and  the  remaining 

tans,  Bfl  to  axe  static.  Define  V^^vnd  V^^^in  the  same 
(analcgcus)  manner.  Ihus, 
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(A-6) 


